Polynomial growth in branching processes 
with diverging reproductive number 
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We study the spreading dynamics on graphs with a power law degree distribution p^ ~ fc ' with 2 < 7 < 
3, as an example of a branching process with diverging reproductive number. We provide evidence that the 
divergence of the second moment of the degree distribution carries as a consequence a qualitative change in 
the growth pattern, deviating from the standard exponential growth. First, the population growth is extensive, 
meaning that the average number of vertices reached by the spreading process becomes of the order of the graph 
size in a time scale that vanishes in the large graph size limit. Second, the temporal evolution is governed by a 
polynomial growth, with a degree determined by the characteristic distance between vertices in the graph. These 
results open a path to further investigation on the dynamics on networks. 
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Branching processes model the evolution of populations 
whose elements reproduce generating new elements fil 0], 
such as a population of physical particles |3|, |4{], cells |2|], or 
infected individuals |;5J]. A key magnitude determining the 
dynamical evolution of the population size is the average re- 
productive number R, giving the number of secondary parti- 
cles generated by a primary particle. When R < 1 the aver- 
age number of new elements decreases exponentially, while 
it grows exponentially when 1 < R < 00 (lj. On the other 
hand, it has been recently found that R may be unbounded 
for branching processes taking place on graphs with a power 
law degree distribution |5l|fj|7j[8l|9|], where by unbounded we 
mean that R diverges with increasing graph size. This obser- 
vation is extremely important since several graphs represent- 
ing interactions among human or computers are characterized 
by a power law degree distribution 11(1 llli Il2l 1 1 3l Il4ll . re- 
quiring us to consider branching processes with an unbounded 
average reproductive number. 

Barthelemy et al fl5ll have recently studied the spreading 
dynamics of an infectious disease on a graph with a power 
law degree distribution. Using a mean-field approach they 
obtained that the average number of infected vertices grows 
exponentially in time with a characteristic time r ~ (fc 2 ) -1 , 
where (k 2 ) is the second moment of the degree distribution 
Pk- For graphs where pk ~ k 1 with 2 < 7 < 3 the second 
moment diverges and r — ► with increasing graph size, pre- 
dicting that all vertices will be instantaneously infected 1 15]. 
A disease that spreads at constant rate, however, cannot spread 
to all vertices in a time scale much smaller than the inverse of 
the spreading rate, indicating that the predicted exponential 
growth should not dominate the system's dynamics. 

In this work we study branching processes with an un- 
bounded average reproductive number using a spreading pro- 
cess on a graph as a case study. When the degree distribution 
has the power law tail pk ~ k^ 1 with 2 < 7 < 3 we ob- 
tain that the exponential regime is followed by a polynomial 
growth in time, a result that is completely unexpected based 
on previous mathematical studies. We also show that both the 



characteristic time separating the exponential and polynomial 
regimes and the polynomial degree depend on the characteris- 
tic distance between vertices. More important, in the limit of 
infinite graph sizes the exponential regime is virtually absent, 
indicating that the polynomial regime is a novel and character- 
istic feature of the spreading dynamics on graphs with degree 
exponent 2 < 7 < 3, and more generally of branching pro- 
cesses with an unbounded average reproductive number. 

Consider a spreading process on a graph with a tree-like 
structure. At t = a vertex selected at random is infected by 
a "virus", which can then propagate to other vertices through 
the graph edges. The causal tree representing the spreading 
process can be modeled as a branching process. Each vertex 
in the causal tree represents an infected vertex in the original 
graph, and each arc in the causal tree represents the generation 
of a secondary infected vertex from a primary infected vertex. 
The out-degree of a vertex in the causal tree gives the num- 
ber of other vertices it infects, i.e. its reproductive number. 
In turn, the length of an arc A^B in the causal tree gives the 
generation time, the time elapsed from the infection of the pri- 
mary case A to the infection of the secondary case B. Finally, 
the vertex generation coincides with the topological distance 
from the first infected vertex, the root, in the original graph. 

We assume that the reproductive numbers are independent 
random variables with the probability distribution q[ an d av ~ 
erage reproductive number R^ = £\ q[ d 'k, parametrized 
by the generation d. The parametrization by d is introduced to 
take into account that the degree distribution may change sig- 
nificantly from generation to generation 00. We also as- 
sume that the generation times are independent random vari- 
ables with the distribution (r) and the probability density 
5 ( d )(r) = dGW(r)/dT. Let P$\t) be the probability dis- 
tribution of the number of vertices N that are found at time t 
in a branch of the causal tree, given that branch is rooted at a 
vertex at generation d. Because of the tree structure we can 
write the recursive relation 
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with the boundary condition Pjy (t) = §n,i> where D is the 
maximum distance between two vertices on the graph. The 
sum over k runs over the possible reproductive numbers of the 
reference vertex, while the sum over Ni, i — 1, . . . , k, runs 
over the possible number of infected vertices in the branch 
rooted at the ith neighbor of the reference vertex. These sums 
are then restricted by the Kronecker delta to configurations 
satisfying 1+Xa=i -^i = N. Finally, within the [• • •] we have 
the probability that the branch rooted at the z-th neighbor has 
Ni infected vertices at time t—T, averaged over the generation 
time distribution G(r). The product structure in Q suggests 
the use of the generating functions 
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for the reproductive number and the number of infected ver- 
tices, respectively. From Q-0 we obtain 
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with the boundary condition F( D > (x) — x. From this equa- 
tion we obtain the average number of infected vertices on the 
branch rooted at a vertex at generation d, 



N^ d \t) 



dF( d \l,t) 
dx 



(5) 



with the boundary condition N^ D \t) = 1. Iterating this 
equation from d = D to d = we obtain the average 
number of infected vertices at time t, N^ Q '{t), and the av- 
erage number of new vertices infected between t and t + dt, 
n(t) = dN^(t)/dt, resulting in 



(6) 



Next we consider the cases when: (/) the reproductive num- 
ber of vertices other than the root has the same statistical prop- 
erties, i.e. = R and = R for d > 1, and (ii) the 
infection is transmitted from an infected vertex to a suscepti- 
ble (not yet infected) vertex at constant rate 1/Tq. This last 
assumption corresponds to an exponential distribution of gen- 
eration times G^ d \r) = 1 — exp(— t/To), with average gen- 
eration time Tq. Under these approximations from (|6j and @ 
we obtain 
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The sum in (|8jl is the Taylor series expansion of exp(Rt/To), 
up to the D—l order. It actually approximates an exponential 
function depending on the ratio of t/ro, where 



where 
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is the average number of vertices at generation d, and the sec- 
ond factor is the probability that the infection has reached a 
vertex at generation d, where * denotes the convolution oper- 
ation, for instance gW * grW (t) = f* drgW (r)gW (t-r). 
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When f < T we obtain 
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becoming an exponential growth for R > 1 OJ Il5ll . In con- 
trast, when f > To we obtain a polynomial growth followed 
by an exponential decay: 
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In general the time scale To depends on the graph size 
No. For random graphs with an arbitrary degree distribution 
< ?i° ) = Pk and q^' = (k-l)pk-i/{k) for d > 01, resulting 
in R ~ (k) and i? ~ (fc 2 ), where (k) and (fc 2 ) are the first 
and second moments of the degree distribution. In this case 
we obtain the following scenarios: 

(0 When the tail of the degree distribution decays faster 
than pk ~ fc~ 3 the diameter scales as D ~ log Ao |7], while 
R is constant or approaches a constant in the large graph size 
limit. Thus, from l|9} it follows that 



to ~ ^-logTVo . 
R 



(12) 



tq, where 



In this case the exponential growth last till t 
To — > oo when ATq — > oo. 

(z'z) When the degree distribution has the power law tail 
Pk ~ fc -7 with 2 < 7 < 3, the diameter D increases at 



most as log A*o while R 

from (|9} it follows that 
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Thus, 



(13) 



The initial exponential growth is thus a finite size effect re- 
stricted to t <C To, where To — > when A*o — » oo. Following 
this vanishing time window the number of infected vertices 
is already of the order of the graph size No {RR ~ A*o) and 
its temporal evolution is polynomial (II H . with a degree de- 
termined by the characteristic distance between vertices in the 
underlying graph. 

To check the validity of our calculations we perform numer- 
ical simulations of the susceptible infected (SI) model on ran- 
dom graphs with a power law degree distribution pf. = Ak _1 . 
Within this model, vertices can be in two states, susceptible or 
infected, and infected vertices transmit the infection to each of 
its neighbors at a constant rate 1/Tq l2oll . We generate ran- 
dom graphs with a power law degree distribution using the al- 
gorithm proposed in 12111 . Then we generated single outbreaks 
on these graphs starting from one infected vertex. Finally, we 
take averages over 10,000 outbreaks starting from randomly 
selected vertices, and over 100 graph realizations. 

When 7 > 3 the spreading dynamics is better described by 
an initial exponential growth (Fig. ^), in agreement with dlOt 
and previous mathematical approaches |15, 20, 22]. In con- 
trast, when 2 < 7 < 3 the spreading dynamics is better de- 
scribed by (Fig- E 5 ), and the exponent D resulting from 
the fit to the numerical data scales linearly with the average 
distance between nodes (see inset of Fig. \ljp). In a more real- 
istic scenario, we use the SI model to simulate the spreading of 
a routing table error on the Autonomous System (AS) network 




FIG. 1: Fraction of infected nodes p(t) — n(t)/No as a function of 
time resulting from SI model simulations on random graphs with a 
power law degree distribution pt = Ak" 1 , with 7 = 3.5 (a) and 
7 = 2.5 (b). Different symbols correspond with different graph 
sizes: iV = 1000 (circles), 10,000 (squares) and 100,000 (trian- 
gles), (a) For 7 = 3.5 the spreading dynamics is characterized by an 
exponential growth (line), as predicted by 1101 . (b) For 7 = 2.5 the 
number of new infections is better described by WW (line). There 
are some deviations at short times, but they get reduced with increas- 
ing the graph size. The inset shows the exponent D resulting from 
the fit of 1 1 1 i as a function of the average distance (d) between two 
nodes in the graph. The increase in (d) is obtained by increasing 
the network size from No = 1000 to 10,000, and 100,000. The line 
emphasizes the linear scaling between D and (d) . 



representation of the Internet |23]. This network is character- 
ized by a power law degree distribution with 7 ~ 2.1 full , 
but it also exhibits degree-degree correlations II ill and a large 
degree dependent clustering coefficient l l24ll . Yet, the average 
number of new infections is well fitted by ( II Q , indicating that 
our predictions are also valid for graphs that are not random 
as well (see Fig.|2j- 

With relevance to the spreading of computer virus and 
worms among email users, there is empirical evidence indi- 
cating that Email networks are characterized by a power law 
degree distribution with 2 < 7 < 3 1 13, The transmis- 
sion rates of computer viruses are, however, of the order of 
their typical detection times, making difficult the empirical 
observation of the initial epidemic growth. With relevance to 
sexually transmitted diseases, there are several reports indi- 
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FIG. 2: Fraction of infected nodes p(t) = n(t)/No as a func- 
tion of time resulting from SI model simulations on AS networks, 
of September 1997 with No = 3015 (circles) and of October 2001 
with No — 10515 (squares). The line is a fit to lilt resulting in 
D = 4.7 ±0.1. 

eating that the network of sexual contacts is characterized by 
a power law degree distribution fl^.E3.E3l . with an exponent 
7 > 3 for some communities and 2 < 7 < 3 for others. This 
fact together with the results obtained in this work represent 
a possible explanation for the observation of both exponential 
and polynomial HIV epidemic growth in different populations 
I22t l27l 12811 . The available data is, however, not sufficient to 
make a definitive conclusion. 

In a more general perspective, our results indicate that the 
degree statistics is not sufficient to characterize the spread- 
ing dynamics, and probably other dynamical processes, tak- 
ing place on graphs with a power law degree distribution with 
exponent 2 < 7 < 3. To determine the characteristic time tq 
and the polynomial degree we need the characteristic distance 
between vertices in the graph as well. The amount of infor- 
mation needed to determine the distance between vertices is, 
however, more difficult to collect, in principle requiring the 
complete mapping of the graph. In this respect, the devel- 
opment of realistic graph models that can accurately represent 
the real graphs will be extremely valuable, allowing us to char- 
acterize the distance statistics from the degree statistics. 
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